Asymptotic freedom and IR freezing in QCD: the role of gluon 

paramagnetism 

Yu. A. Simonov 1 '!*! 

Institute of Theoretical and Experimental Physics 
Moscow, Russia 

■ Paramagnetism of gluons is shown to play the basic role in establishing main 

properties of QCD: IR freezing and asymptotic freedom (AF). Starting with Polyakov 

o 

■^j- ■ calculated and shown to ensure not only the classical result of AF but also IR freezing. 



X 



For the latter only the confining property of the background is needed, and the 
effective mass entering the IR freezing logarithms is calculated in good agreement 



Q^' with phenomenology and lattice data. 

^ i 1. INTRODUCTION 

OO 
CO 

lO ' The notion of asymptotic Freedom (AF) is basic in establishing QCD as a selfconsistent the- 
ory [1]. The extrapolation of the QCD coupling constant a s (Q) to larger distances (smaller 
momenta Q) leads however to inconsistencies of several kinds in the pure (nonbackground) 
perturbation theory: 



1. The appearance of Landau ghost pole (and other singularities in higher orders) pre- 
cludes extrapolation to small Q 

2. IR renormalons make the whole perturbation series not summable even in the Borel 
sense 



Dation series in the Minkowski space-time has difficulties and 
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3. The treatment of pertur 
should be reformulated 

At the same time the IR behavior of a s in experiment |5j and on the lattice 6] does not 
show irregularities in the Euclidean region, Q 2 > 0, and is compatible with IR freezing. 



*Electronic address: 



simonov@itep.ru 



To ensure this nonsingular behavior a special type of theory was suggested 



0. 



iminating 



Landau ghost pole from the beginning, which is phenomenologically successful |8]. 

To understand what dynamical mechanism makes QCD perturbation theory consistent 
and brings in the IR freezing, as seen on lattice and experiment, the Background Pertur- 
bation Theory (BPTh), formulated earlier in (9], was considered, treating background as a 



d with the property of confinement 10]. 



11], that the basic effect of this confining background is to make a s (Q 2 



strong collective fie 
It was shown in 

finite at all Q 2 > and thus precluding appearance of Landau ghost pole and IR renor- 
malons. Moreover, the AF logarithm approximately keeps its form at large Q 2 , In — , but 

QCD 

at small Q 2 the argument acquires the additional term, which looks like the two-gluon mass 

nn 

. This type of form was suggested before [12J, H3j, however in QCD 



Mf yielding In 



Q 2 +Ml g 



"qcd 



the appearance of gluon mass is forbidden by gauge invariance. As will be seen the term 
M 2g actually has the meaning of the two-gluon mass, where gluons are connected by the 



adjoint string. In 



11| the IR freezing was considered in the framework of the static QQ po- 



tential, and the exact general form of IR behavior and exact value of Mi g were not actually 
given. In the present paper we present a more general derivation of the IR freezing based 



on the Polyakov background approach 14j, where the basic one-loop element is the scalar 
self-energy (gluon loop) operator U(Q 2 ). As will be shown, in the confining background 
H{Q 2 ) acquires the two-gluon mass M 2g ~ 2 GeV, and ensures both AF at large Q 2 and 
IR freezing at small Q 2 . The explicit value of this mass is estimated and it is shown, that 
numerically IR freezing is not universal: the IR freezing behavior and mass (to be called 
IR mass) depends on the embedding process. For comparison the background perturba- 
tion theory for the static QQ system is considered in the one-loop approximation and it is 
shown, that the corresponding IR mass is much lower, M 2 ? ~ 1 GeV 



his latter value is 



in good agreement with phenomenological description of IR freezing 15|, [l6| as well as with 
lattice determinations of a s Q]- The generalization can be considered as well, leading 
to the inclusion of multigluon states in the asymptotics of a corresponding Green's function, 
which coupled by confinement. Thus all theory becomes finite and devoid of IR renormalons 
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191 ]. hence well defined in the Euclidean region. 



To go beyond Euclidean region, one needs to define better the singularity structure of 
the perturbative series. The logarithmic singularities of the free PTh are not physical as 
well as those in BPTh. To simplify matter one can take the limit iV c — > oo, where all QCD 



3 



amplitudes contain only poles [20 ]. 

The corresponding extrapolation was done in 



21 



22] where it was shown, that equidistant 



mass squared spectra of hadrons allow to replace all logs by Euler -^-functions, and thus 
obtain for finite Q 2 < simple poles in Minkovskian region, while for large Q 2 > in 
Euclidean region one has standard logarithmic terms. The whole scheme works nicely for 
both (3(a s ) and a s and agrees well both with lattice and phenomenology 2l|, |22 ] . 

In all these considerations the nonpositive definiteness of the /9-function of SU (N c ) theory 
is crucial, and the latter is due to gluon paramagnetic terms in Lagrangian and gluon Green's 
function. The plan of the paper is as follows. Section 2 is devoted to the extrapolation of 
the Polyakov method to the IR region. Section 3 contains similar treatment for the static 
potential QQ system. In section 4 summary and discussion of results is given. 



2. ONE-LOOP EVOLUTION OF a s BY THE POLYAKOV METHOD 

As in [jjj], one starts with the gluonic action S = J F^F^ u d 4 x, defined at the scale R\ 
(momentum scale Ai = 1/Ri) and consider Wilson transformation to the scale i?2(A2), which 
can be considered as the change of the effective integral volume in S from Rf to R\. Sepa- 
rating gluon field into valence gluons a M and backgound B^. 

= a M + Bp, (1) 

one can expand in a M , keeping quadratic in terms 

« = < ((DlUal - 2gF« u (B)a b IM atr bc ) . (2) 

As Polyakov mentions, the first term, proportional to D 2 , gives rise to diamagnetic inter- 
action of valence gluon with background, ga'j 1 B , ld\alf abc , while the second term is param- 
agnetic interaction of gluon spins with background. Both can be expressed in second order 
through the scalar gluon self-energy II (x — y), which corresponds to the loop diagram of two 
massless scalars, and in case of no background is 

U (x) = G 2 (x) = ^-L-j. (3) 

The resulting expression for the change of one-loop correction from the scale R± to the scale 
i?2 is 

1 1 +^f(Ri, R 2 )X = ^N C (4) 



g 2 (R 2 ) g^R,) 4tt' v " * n u 3 
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where we have defined 

"\x-y\=R2 



f(R 1 ,R 2 )= I ' d\x-y)Yi{x-y) (5) 

J \x—y\=R\ 

which yields the standard expression in the free case (no background), 

/ ( J R 1 ,i? 2 ) = -^ln^. (6) 



It is this behavior, which produces AF at small Ri and Landau ghost pole |2[ appears when 
the r.h.s. of (jlj) vanishes. In case of nonzero background it is necessary to take into account, 
that "scalar gluon" propagator in background, G(x) is no nore free and massless. Moreover, 
if one takes into account confinement, then the product G 2 (x) should be replaced by the 
two-gluon white Green's function, i.e. the two-gluon glueball Green's function. G 2g (x), and 
the resulting evolution function 

f(R 1 ,R 2 ) — > f 2 g(Rl, R2), 

/■\x\=R2 
d'xU^x). (7) 
,x\=Ri 

It is our purpose below in this section to calculate f 2g both in coordinate and in the 
momentum space, proving the IR freezing in an explicit way. 



To proceed we shall use the exact Fock-Feynman-Schwinger Representation (FFSR) 23], 
for the 2g Green's function n(x, y) in the nonzero background 

/•oo /*oo 

U(x,y) = j Q d Sl j ds 2 (Dz^) xy (Dz^) xy e- K '- K2 W a (x,y) (8) 



(0 



4 JO \ dt 

insertions, 



where Ki = | J Si ( ) dr^ and W a is the Wilson loop with paramagnetic gluon spin 



W a (x,y) = Pexp \ ig j^^ApdzJj exp (tog jf* F(z(r))^ dr. (9) 

Here C (x, y) is the loop contour formed by the paths of two gluons from the point x to the 
point y. Averaging over the vacuum configurations one obtains U.(x,y), which is_ expressed 
only in terms of einbein parameters to be found from the 2g Hamiltonian 
Appendix 3 of 



24 



25|, (see 



11| for details) 



M1+M2 rp 
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where we have defined T = \x — y\,jl 



G(0,0,T) = (0|e-^ r |0) = El^£ ) (0)re 



1 2 -Ei n) T 

2g 



(11) 



n=0 



Here B.2 g is the 2g Hamiltonian with confinement and spin-dependent interaction, derived 



in E24, 



25| . and "0" refers to zero intergluon distance in the initial and final state. 



The S'-wave spectrum of the Hamiltonian to the lowest order in spin splittings is well 
known 24-2f| 

jx 2/3 



M<^> = f*±*i + e n Qi), e n {H) = (2~^a^a(n), a(n) ~ 



n H — 
2 



(12) 



Here cx a dj = |o" found is the gluonic string tension, and it is conceivable, that this gluonic 
string does not decay for N c — > 00, however even for finite N c the main results are not 
sensitive to the high excitations, as will be seen, and hence to the string decay. Inserting 



(n) 



E 2g (n) = M n from PJ and 



4?r ' 



[271 ] . one obtains 



(13) 



where A 



4(5hrW/ 2 



Following [rjj, for large T one can do integration over dfXi, d\i 2 using the steepest descent 



method, which yields stationary point fi^ = ^ = jM n , where 



a(n) 



3/4 



(14) 



and the resulting form of Tl(x, y) is 



T 3 / 2 



-M n T 



(15) 



At small T one instead goes from the sum over n to the integral, which yields 



U AF (x,y) 



1 



_Ml+M2 T 



16vr 4 T 4 



(16) 



4(2vr) 4 T 2 Jo 

which reproduces the free result ([3]). 

Let us now turn to the momentum space. The analysis of Polyakov in the free case can 
be written in the form 

1 1 



9 2 (Q 2 ) 9 2 {ti) 



-n para (Q 2 ) + n dia (Q 2 ), 



(17) 
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where 



n para (Q 2 ) = 4iv c n(Q 2 ), n dia (Q 2 ) = ^n(Q 2 ), (is) 



and I1(Q 2 ) in the free case is simply a scalar gluon loop, which after renormalization takes 
the form 

IWQ 2 ) = / j= , 4 tr-nv = ~ih ln (%) ■ ( 19 ) 

Jfio (27r)V(P + Q) 167r \ A*o / 
Let us now turn to the case of perturbation theory in the confining vacuum. In this case 

two gluons in the loop form bound states, and we can use the spectrum, given in (|14p for 

large n 

lTia a ( \\ 2 f 2 \ 2 4vra a 



M » = vr t" + = m + c = 47R7 ° U J • m "vr (20) 



For the WKB spectrum in the linear potential a a r one would obtain instead 



28) 



m WKB = 27rcJ a, C W KB = 47TCT a (21) 

and we shall exploit these values (differing by 15% from m 2 and c respectively) in what 
follows. 



In terms of the bound glueball states H(Q 2 ) — » n conf ((5 2 ) can be written as [29J] 



= yj = i ^> ( o)|» = £ (22) 

Replacing in ( l2"2"j) the sum over n by the integral and renormalizing the integral in the same 
way as in ( fl9l) one obtains 



n conf (Q 2 , „ 2 ) = --^ In g2 + ^ KB (23) 
1071 y^Q 



This latter form coincides at large Q with the perturbation theory result 

n conf (Q 2 ,/i 2 ) n , =nUQ 2 ) = -T^- 2 ^(%) (24) 

Thus the charge evolution to the leading order in the confined vacuum can be written as 

Tm^TW)- 1 ^^ 2 '^ (25) 

One can see from (|2"7j) . that for large Q ^~ 2 m , n conf (Q 2 , /i 2 ) < and the AF appears, 

s 2 (Q 2 ) < <rV)- 
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For nonasymptotically large Q 2 , but still when Q + 2 m is large enough, the standard form 
of one- loop result for a s {Q 2 ) using fl23|) . fl25l) acquires the form 

^(Q 2 ) = (26) 

This form coincides with the one, proposed long ago in [3], where m was associated 
with the effective mass of two gluons, As one can see, this notion of mass is in reality 
extended to the ground state mass of two gluons connected by the adjoint string, i.e. a 
ground-state glueball mass M 2fl (0 ++ ). 

Numerically, however, M 2g (0 ++ ) is large, from (12TT) M 2g (0 ++ ) = iti^kb = 1-6 GeV, which 



agrees with explicit calculations in 



24 



251 ]. while phenomenological estimate for m in the 



a s , entering the static QQ potential, is m ~ 1 GeV 15|, ll6j. In the next section we consider 



this situation in detail and shall find rngQ for the static potential. 

3. ONE-LOOP EVOLUTION OF a s FOR THE STATIC QQ POTENTIAL 



The purely perturbative derivation of static potential is given in 30]; for the case of con 



finement this situation was considered in detail in [HI]. Below we shall give the main results 



of [Uj and, as a new element, we estimate numerically the IR freezing mass m in Eq. ( 126]) 
for the explicit case of the static QQ potential, to be called ttlqq. 

One starts with the Wilson loop with rectangular contour C of size R x T containing 
both nonperturbative confining background and valence gluons a M . Expanding in powers 
of (ga^), one obtains a series of diagrams with valence gluon exchanges in the background 
field and after vacuum averaging 

(W(B + a)) B , a = exp(-V(R)T - perimeter) (27) 

one has terms 

V(R) = V q (R) + a s V 2 {R) + a 2 s V A {R) + ... . (28) 

Here V q (R) = afR at large R, while V n (R),n > 2, contains up to n gluons propagating 
inside the minimal surface S bounded by the contour C . 

Thus Vi{R) corresponds to the one-gluon exchange; while V±{R) contains a gluon loop 
(minus ghost loop) on the gluon propagator, a triangle vertex part and double gluon exchange 
and we take the limit N c — > oo, so that each gluon line is represented as a double fundamental 
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line. We also take the limit T — > oo to define static potential properly. In this case the 
change in the area of the minimal surface 5* due to gluon propagation is only due to inner 
closed loops in V A (R), while V^-R) and all single gluon lines are unaffected by confinement, 

dx A dy A a^C 2 (f)T 



T/fD\ „2 g 2(/) 



We can write the contribution of V 2 + V4 in the form 



V 2 (R) + V A (R) 



Jo (a?4 - yi) 2 + R 2 

l + a^f(R)), 



R 



R 



where f(R) was computed in case of no confinement in JjJ, [30] 



60 



47T 



R 



f (R) = -±ln - , 5-1///. 



In the MS 1 scheme fo(R) was found to be 31 1 



fo MS (R) = ^{^ 2 R 2 ) + 2i E ) + l 



7T V12 



(29) 



(30) 



(31) 



(32) 



When confinement is included in the background, / con f(-R) is expressed through the gluon 
selfenergy term Il(x — y), introduced in the previous section. 



^conf(^) 



fconf(R) 



d 4 rll(r) 
|R- r| 



7 11 Ar 

&n = — N r . 



R An 2 

In the momentum space the Fourier transform of the / CO nf / R can be written as 

t (Q)=/<i^„, t (fl)e-Q R =^S2), 

J 7T U 



Vconii 

and using (1231) this can be written as 

^conf(Q) 



i67r 3 cr 



111 



Q 2 + m 2 
/i 2 



Hence the total one-loop potential in momentum space has the form 



V 2 (Q) + Vt(Q) 



C 2 (f)a 



(o) 



1 



b , _ Q 2 + 



47T 



-a. In 



(33) 



(34) 



(35) 



(36) 



4tt 2 Q 2 

It is now essential, that the string, connecting the gluons in the internal loop, is fundamental, 
and therefore 



m 



QQ 



27T0-/, m Q Q = 1.06 GeV. 



(37) 



This is important result, since the IR freezing of the gluon-exchan ge p otential is essential 



15 



16 



32|). as well as in 



for the quark model calculations, e.g. of hadron masses (see 
different QCD processes; it can also be tested on the lattice. In the next section we shall 
compare the result of fl37|) with other approaches. 
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4. RESULTS AND DISCUSSION 

We have studied in previous sections the a s renormalization with and without confinement in 
two different settings: the Polyakov background setting in section 2 and the QQ interaction 
in section 3. We have found, that in the quenched (N c — > oo) case the one-loop a s is given 
by the same equation in all cases 

An - 11 

*s(Q 2 ) = = Q2+m2 , b = -N c . (38) 
Oo m — j2 — ° 

Here m 2 = for the case of no confinement in both types of setting, confirming that the AF is 
a universal phenomenon. However, the IR mass m is not universal, it is m 2 = m 2 g = 2na a in 
case of Polyakov background approach and m 2 = rn 2 Q Q = 2iraf in case of the QQ potential, 
and the latter can be written as 



V QQ (Q) = -AnC 2 (f)^±, (39) 



and ay(Q) to one loop is the same, as in ( 13 8 p with A — > Ay, while in the two-loop approxi- 
mation can be written as 1^, Q 

with t B = In h = 102, b Q = 11. 

" Ay 1 7 v 

The coordinate-space representation Vqq{t) was studied in detail in [18j and it was shown, 
that to a reasonable accuracy (better than 10% for r > 0.2 fm) ay(r) can be approximated 
by ay(Q = 1/r), so that 

C 2 (f) 4t _ C 2 (/)_ (1) , 

In the two-loop case Vqq(t) was computed in [l^, Q] and is given by the same Eq.(j4"Tl) 
where now the two-loop a v (r) is [18| 

4 2> (r) = 4"(r) jl + ftW^ + ft(r) (^>) j , (42) 

where JBi(r) = d + 27i(r)6 , 5 2 (r) = a 2 + 271 (r) (61 + 26 ai) + &o72( r ) an d 

2 /- 00 sinx ~ ~ / l + m 2 r 2 \ 

7n(r) = - / dx (t(x)) n } t(x) = In 



7r jo x \ x 2 + m 2 r 



2^2 / ' 
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31 



«2 



4343 



7T 



22 



+ 4tt z - — + — C(3) 
162 4 3 SV ; 



One can compare the one-loop results in p space 



or two-loop result (|4"U|) . and the 



corresponding x-space expressions (l4"Tj) and 



for ay(r) from 



33] were compared with aty (r) i n |l7] and shown to be in good agreement 



with lattice and experiment. Lattice data 



rr 



in the region 0.04 fm < r < 0.4 fm, measured in [33], while the purely perturbative «p er t(^) 
strongly deviates from lattice data already for r > 0.08 fm. a similar good agreement can be 
deduced, comparing a v (r) to the Schroedinger functional lattice method (second reference 
in j^]). On the phenomenological side, hadron spectra and especially fine-structure splittings 
are sensitive to the behavior of acy(r). The saturated (frozen) value of ay = a CT n at large r 



was assumed in the detailed calculations in 
in 



17 



32] and this value of a cr i t agrees well with found 



18] . Moreover, the analysis of the splittings between low-lying levels in bottomonium 



151 ]. yields a cr ; t = 0.58 ± 0.02 in striking agreement with [32] . 

A detailed analysis of bottomonium splittings in comparison with lattice data and ex- 
periment [k| proves, that the potential Vqq(t) ( 14T1) with ay\r) obtained in the confining 
background perturbation theory with mqq = 1 GeV is in good agreement with experiment. 
This confirms the agreement between the calculations of the present paper and the physical 
reality. 

Concluding one should stress the crucial role of gluon paramagnetism in creating the prop- 
erties of AF and IR freezing. The correct sign of the logarithmic term in ( 126]) is important for 
its Minkowskian extrapolation in the form of the sum of pole terms in the ^-function, which 
in its turn correspond to correct physical poles in a s , as shown in 2l|. Moreover, as shown 



in 



34] , the gluon paramagnetism is responsible for the correct physical behavior of the field 



correlator (F° {x)F® (x)), implying nonzero and positive value of gluonic condensate. 

The author is grateful to A.M.Badalian for useful discussions; financial support of the 
RFBR grant 09-02-00629a is acknowledged. 
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